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ABSTRACT 


A  theoretical  model  was  developed  for  antigen-antibody 
interaction  in  a  receptor-carrier  based  system.  Many 
models  have  been  developed,  both  kinetic  and  equilibrium, 
for  receptor-to-ligand  interaction  in  the  immune 
response.  The  immune  response  is  extremely  complex,  with 
very  large  numbers  of  receptor  sites  on  many  distinct 
cell  types  in  the  system.  Therefore  the  system  needs  to 
have  many  restrictions  upon  it  to  make  it  manageable. 

Most  models  assume  the  area  for  interaction  is  an 
infinite  sheet  (i.e.  an  endless  two  dimensional  cell 
membrane) .  These  models  are  limited  to  only  discussing 
the  average  number  and  types  of  interactions.  This  model 
is  receptor-carrier  based  which  counts  interactions  on  a 
cell/vesicle  surface  to  allow  questions  about  actions  and 
conditions  per  cell  to  be  addressed.  The  models  which 
have  been  addressed  are  the  ••sandwich"  model  which  takes 
a  two-step  interaction  with  a  linking  ligand  and  then  the 
antibody,  and  the  agglutination  model  which  addresses  the 
inter-cellular  binding  between  two  adjacent  cells  through 
a  bridging  ligand.  Results  show  that  the  sandwich  model 
and  the  agglutination  model  are  less  sensitive 
immunoassay  techniques  than  the  direct  binding  assay, 
although  these  models  are  very  useful  for  describing  more 
complex  interactions. 
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I .  INTRODUCTION 


The  immune  response  in  humans  is  a  highly  organized 
and  very  specific  defense  system  [1].  The  immune  system 
defends  the  body  against  foreign  molecules,  cells,  or 
viruses,  which  are  called  antigens.  The  human  system  not 
only  recognizes  non-human  cells,  but  also  those  from 
other  individuals,  which  explains  the  difficulty  in 
transplanting  organs  from  one  person  to  another  [2]. 

There  are  two  main  types  of  cells,  lymphocytes, 
which  are  involved  in  the  immune  response:  T-cells 
(which  mature  in  the  thymus)  and  B-cells  (which  mature  in 
bone  marrow) .  The  humoral  immune  response  is  produced  by 
the  B-cells  and  is  most  effective  against  pathogenic 
(disease  causing)  cells,  viruses  and  free  toxins.  The 
cell-mediated  immune  response  uses  the  T-cells  and  is 
effective  against  larger  cells,  such  as  parasites, 
cancer,  and  infected  cells. 

The  cell-mediated  response  (Figure  1)  is  more 
complex  than  the  humoral  immune  response  (Figure  2) . 

This  is  the  system  which  is  attacked  by  the  AIDS  virus. 

In  many  responses  the  antigen  is  picked  up  by  a 
macrophage  and  is  presented  to  a  "virgin”  T-cell  which 
initiates  a  response.  The  simplest  T-cell  response  would 
then  produce  the  cytotoxic  T-cell  which  kills  the  antigen 
bearing  cell  and  then  looks  for  other  antigens.  The 


Humoral  Immune  Response 


Phagocytosis 
by  macrophages 

Figure  2 


other  response  initiated  by  the  T-cell  involves  helper 
T-cells  and  suppressor  T-cells.  The  helper  cells  secrete 
lymphokines  and  interleukins  which  activate  macrophages 
and  mast  cells  in  the  area  to  ingest  the  antigen.  The 
helper  cells  can  also  induce  local  plasma  cells  to 
produce  antibodies.  The  helper  cells  also  tell  the 
antigen  bound  T-cells  to  reproduce.  At  the  same  time 
suppressor  T-cells  are  doing  the  opposite  of  the  helper 
cells  in  a  dynamic  balance  to  induce  an  appropriate  level 
of  response.  An  appropriate  response  is  necessary  since 
too  little  response  will  not  rid  the  body  of  the  antigen 
and  over  response  is  what  causes  side  effects  like 
sneezing,  coughing  and  hives. 

The  B-cell  response  is  more  direct  and  simpler  to 
understand.  The  B-cells  move  freely  through  the  blood 
and  lymph  system  of  the  body.  When  a  B-cell  encounters  a 
specific  antigen,  it  will  begin  to  grow,  and  divide. 

This  division  produces  two  types  of  cells,  the  memory 
cells  and  the  plasma  cells.  The  memory  cells  allow  a 
more  rapid  response  to  be  initiated  if  the  same  antigen 
were  to  attack  again.  The  plasma  cells  actually  produce 
antibodies  which  are  used  to  attack  the  antigen.  Figure 
3  shows  the  layout  of  a  typical  IgG  antibody.  There  are 
more  kinds  of  antibodies  than  just  IgG,  but  IgG  is  the 
most  common,  with  serum  levels  up  to  12.0  mg/ml.  The 
antibodies  are  made  up  of  two  light  and  two  heavy  protein 


Antibody  Structure 


Figure  3 


chains.  Antibodies  are  very  specific  in  their  binding 
sites  (Fab  portion) ,  and  hence  differ  in  the  amino  acid 
sequencing  in  that  region,  but  all  contain  an  Fc  tail 
which  is  constant  in  amino  acid  sequencing  from  one 
antibody  molecule  to  another.  This  Fc  tail  is  the 
portion  of  the  antibody  which  helps  initiate  the 
complement  system.  This  system  is  shown  in  Figure  4. 

The  chain  of  proteins  forms  a  pathway  to  destroy  the 
antigen,  and  acts  as  a  catalyst  to  promote  larger 
response  further  along  in  the  process.  The  system  has 
two  methods  of  ridding  the  body  of  the  antigen: 
macrophages  in  the  area  can  be  activated  to  engulf  the 
antigen  or  to  form  a  hole  in  the  cell  membrane  of  the 
antigen  to  lyse  the  cell.  The  humoral  response  can  also 
produce  K  lyphocytes  which  mark  the  antigen  and  kill  it 
through  an  unknown  mechanism  [2]. 

The  purpose  of  this  paper  is  to  describe  in  detail 
the  derivation  of  a  theoretical  model  for  cellular 
interaction.  This  interaction  occurs  in  many  places  in 
the  immune  response;  for  example:  antibody  to  antigen, 
macrophage  to  antigen,  and  helper  T-cell  to  plasma  cell. 

Actual  modeling  of  the  B-cell  response,  considering 
the  huge  number  of  receptor  sites  on  the  B-cell  [3],  on 
an  antigenic  bacteria,  or  on  a  virus  is  unreasonable. 

The  system  needs  to  be  easier  to  control  and  have  fewer 
parameters.  A  more  tractable  system  would  be  to  model  an 
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artificial  system.  The  purpose  of  this  paper  is  to 
describe  in  detail  the  derivation  of  a  theoretical  model 
for  interactions  within  the  immune  system,  its 
usefulness,  and  the  results  of  two  separate  applications 
of  the  theoretical  model. 

The  system  which  is  being  modeled  uses  vesicles  in 
immunolysis  assays,  shown  in  Figure  5  [4].  These  assays 
use  covalently  bound  antigens  on  the  surface  of 
phospholipid  vesicles  [5].  The  vesicles  can  be  made  with 
specific  types  and  numbers  of  receptor  sites.  The 
vesicles  can  encapsulate  markers  like  spin-labelled 
compounds  [6]  or  fluorescent  compounds  like 
carboxyfluorescein  [7].  Carboxy fluorescein  does  not 
fluoresce  in  high  concentrations  like  that  which  is 
encapsulated  in  the  vesicles,  but  when  the  vesicle  is 
lysed  through  the  complement  process  or  by  another  means, 
the  material  is  diluted  and  fluoresces.  This  system  can 
be  modeled  and  a  series  of  questions  can  be  answered  such 
as  the  number  of  antigenic  sites  required  to  activate  a 
response,  and  a  theoretical  limit  to  detection  of  antigen 
or  antibody  with  this  system. 

Another  more  modern  and  exciting  application  of  such 
modeling  is  for  the  system  recently  developed  involving 
immunoadhesins  for  AIDS  therapy  (Figure  6)  [8].  It  has 

been  shown  that  the  AIDS  virus  kills  by  depleting  the 
body's  T4  lymphocytes  (a  type  of  T-cell)  and  the  victim 
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AIDS  Virus  Immunoadhesins 


becomes  susceptible  to  infections  and  neoplasms  (new 
growth,  tumors) .  The  T4  lymphocytes  have  numerous  CD4 
protein  sites  as  receptors.  The  Human  Immunodeficiency 
Virus  (HIV)  envelope  glycoprotein  (gpl20)  has  been  shown 
to  be  strongly  attracted  to  the  CD4  protein.  Most 
therapies  are  ineffective  against  the  AIDS  virus  because 
its  variable  receptors  keep  changing  and  modifying.  The 
one  common  denominator  is  its  attraction  to  the  CD4  site 
on  the  T-cells.  After  the  virus  has  attached  itself  to 
the  T-cell  membrane,  it  is  internalized,  uncoated,  and 
alters  the  cellular  DNA.  It  would  seem  that  an  increase 
in  free  CD4  protein  in  the  blood  would  tie  up  the  coat 
protein  on  the  virus  (gpl20)  and  not  allow  it  to  attach 
to  the  T-cells.  The  problem  with  this  approach  is  that 
CD4  only  has  an  effective  lifespan  of  twenty  to  thirty 
minutes  in  the  blood.  Researchers  have  found  a  way  to 
develop  a  special  antibody,  called  an  immunoadhesin, 
which  has  CD4  receptor  sites  and  an  Fc  tail.  Experiments 
predict  a  lifespan  approaching  21  days,  the  same  as  for  a 
natural  IgG.  The  goal  of  this  program  is  to  tie  up  the 
virus  so  it  cannot  attack  the  T4  cells  and  allow  the 
normal  immune  response,  involving  the  Fc  portion  of  the 
immunoadhesin  and  the  complement  system,  to  kill  the 
virus.  This  system  can  be  represented  by  monovalent 
receptors  (the  gpl20  coat  protein  sites)  on  the  viral 
surface  and  bivalent  ligands,  the  immunoadhesins.  This 


model  can  describe  the  specific  actions  on  the  viral 
surface. 

Many  questions  still  remain  about  the  specific 
methods  and  limiting  characteristics  of  the  immune 
system.  A  theoretical  model  would  be  useful  to  answer 
some  of  these  questions.  Since  the  antibody  to  antigen 
binding  constant  is  ambiguous  [9]  and  a  crucial  factor  in 
explaining  cell  to  cell  interactions,  it  is  a  central 
goal  of  any  theoretical  model  to  find  a  realistic  value 
for  this  constant.  Varying  this  parameter  while  holding 
the  other  conditions  to  those  of  an  actual  experiment 
should  allow  this  constant  to  be  approximated. 

Attachment  to  the  cell  by  first  one  protein  followed  by 
another  is  a  possibility;  modeling  this  two  step  process 
could  show  interesting  effects.  It  is  also  believed  that 
agglutination  of  the  antigens  into  clumps  helps  the 
killer  cells,  macrophages,  and  the  complement  system  to 
work  more  efficiently  [1].  A  model  of  these  situations 
would  help  find  the  system's  limits,  characteristic 
constants,  ideal  conditions  and  number  of  bound 
antibodies  needed  to  produce  cell  lysis. 

The  models  in  this  paper  are  receptor-carrier  based 
(10].  This  form  of  modeling  is  significant  since  it 
allows  questions  to  be  asked  concerning  the  number  of 
antibodies  bound  or  cross-links  per  cell/vesicle.  Other 
models  have  been  proposed  [11,12]  which  are  kinetic  or 


equilibrium  studies  of  monovalent  or  multivalent 
receptors  or  ligands.  These  models  assume  the  surface  of 
interaction  to  be  infinite  so  only  questions  of  types  or 
average  number  of  structures  can  be  asked.  The  counting 
done  by  this  model  keeps  track  of  structures  on  the 
vesicle  surface.  It  is  assumed  that  the  initial  antibody 
attachment  to  the  vesicles  is  caused  by  a  receptor-ligand 
interaction.  This  initial  attachment  rate  is  considered 
to  be  independent  of  number  of  sites  occupied  if  there 
are  sites  remaining  (steric  considerations  are  not  taken 
into  account) .  The  subsequent  attachment  of  multivalent 
antibodies  (in  the  agglutination  model,  the  bivalent 
receptor  model,  and  the  bivalent  ligand  sandwich  model) 
do  take  steric  effects  into  account.  Finally  it  is 
assumed  that  the  receptor  sites  are  able  to  migrate 
freely  [13]  to  allow  adequate  mixing  for  the  equilibrium 
equations  to  be  relevant. 
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II.  GENERIC  MODEL  DERIVATION 

The  models  are  defined  by  a  system  of  conservation 
equations.  The  derivations  are  not  intuitive,  so  an 
example  will  be  provided.  The  simplest  case  is  one  of 
monovalent  receptors  on  the  cell  and  monovalent  ligands 
[14].  This  system  is  represented  in  Figure  7. 

This  system  has  a  vesicle,  or  cell,  covered  by  n 
monovalent  antigenic  receptors.  The  binding  ligands  are 
monovalent,  like  Fab  fragments.  The  system  has  a 
degeneracy,  (i.e.,  number  of  arrangements)  which  is 
defined  as  d^.  Degeneracy  is  the  number  of  equivalent 
arrangements  of  bound  ligands. 

r  n  ^ 

di  =  n!/[i! (n-i) ! ]  =  |  |  (1) 

»>  i  J 

This  is  derived  by  taking  the  number  of  possible 
arrangements  of  receptors  (n!)  and  dividing  by  the  number 
of  arrangements  of  the  different  sites  which  are  on  the 
cell,  including  (i!)  the  number  of  identical  sites  which 
are  bound  and  (n-i) !  unbound  sites.  The  ligands  need  to 
be  specific  for  the  receptor  sites  and  are  governed  by  the 
binding  and  dissociative  rate  constants,  kj^  and  k_^. 

These  constants  are  related  to  K,  the  equilibrium  binding 
constant,  by 


K  =  kj/k.-L 


(2) 
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Monovalent  Receptor-Monovalent  Ligand  Model 


Figure  7 
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The  Law  of  Mass  Action  representing  the  i  th  antibody 
binding  is  represented  by 


K 

vi-l'  +  a  »  Vit  (3) 

which  can  be  represented  by  the  following  arrangement  where  the 
'  represents  one  of  the  specific  arrangements  among  the  total 
d^/  possible  arrangements  of  V-^,  and  "a”  represents  the  free 
antibody  concentration: 

Vi/  -  (aK) Vi_i/  (4) 

The  total  concentration  of  Vi  can  be  defined  as 

Vi/  =  V^di  (5) 

Repeating  equation  (4)  i  times  gives  the  general  equation 

Vi  =  di(aK)iv0  (6) 

where  VQ  is  the  concentration  of  cells  without  any 
antibodies  bound. 

The  total  number  of  vesicles  and  antibodies  in  the 
system  cannot  change.  This  allows  conservation  equations 
to  be  found.  The  conservation  equations  for  the 
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system  can  be  derived  as 

n 

V°  =  l  di(aK)iv0  (7) 

i=0 

n 

Ab °  •  l  id± (aK) XV0  (8) 

i=0 

where  V°  is  the  initial  concentration  of  vesicles  and  Ab° 
is  the  initial  concentration  of  antibodies. 

The  mathematical  identity  of 

n  n 

l  (.)  X1  =  (1+X)n  (9) 

i=0  1 

allows  equation  7  to  be  reduced  to 

V°  =  VD (1+aK) n  (10) 

and  the  derivative  of  equation  9  allows  equation  8  to  be 
reduced  to 

Ab°  =  a  +  nV°aK/ (1+aK)  (11) 


A.  General  Derivation  Summary 


This  is  the  format  which  will  be  used  to  derive  all 
future  models: 

1.  Identify  possible  structures  and  parameters 
governing  receptor  to  ligand  binding. 

2.  Calculate  degeneracy  through  possible 
arrangements ,  including  unbound  sites. 

3.  Strip  (stepwise)  possible  arrangements  to  V0 
(entirely  empty)  arrangement. 

4.  Derive  conservation  equations  for  vesicles  and 
antibodies. 

5.  Address  specific  questions. 


III.  SANDWICH  MODEL 


The  first  model  which  will  be  discussed  in  depth  is 
the  "sandwich"  model  [4].  This  model  is  shown  in  Figure 
8.  This  model  is  different  from  the  monovalent  receptor/ 
monovalent  ligand  model  by  the  ligand  being  bivalent,  with 
one  end  specific  to  the  receptors  on  the  vesicle,  and  the 
other  specific  to  the  monvalent  antibody.  Now  the  vesicle 
receptors  can  have  three  possible  arrangements:  empty, 
with  the  bivalent  "link"  attached  (f),  and  having  the  link 
attached  to  an  antibody  and  the  receptor  site  (s) . 


Sandwich  Model 

Monovalent  receptors-  specific  links 


Figure  8 


The  degeneracy  for  this  model  is  derived  as  in  the 
monovalent  model  and  it  results  in 

df,s  =  n!/[f!s! (n-f-s) !]  (12) 

The  Law  of  Mass  Action  for  the  model  is  represented  by 
K 

Vf,s'  «  vf+l,s-l'  (13> 

where  K  is  represented  by  Equation  2  for  the  s (second) 
step.  By  repeating  this  equation  s  times  the  general 
equation  for  the  concentration  with  f  links  and  0 
sandwiches  can  be  shown  to  be 

vf,s/df,s  =  (^sa) Svf+s,o/df+s,o  (14) 

where  a  is  the  free  antibody  concentration  and  Ks  is  the 
equilibrium  constant  for  that  reaction.  The  Law  of  Mass 
Action  equation 

^f+s, o/df+s , o  =  (Kf4 ) (Ksa) sVf+s-l t o/df+s-l , o 

(15) 

can  be  repeated  f+s  times  to  yield  the  general  equation 


vf,s  -  <Jf,s<Kfl>t+s<Ksa)Sv0>0 


(16) 


is  the 


where  i  is  the  free  link  concentration,  and  VQ  0 
concentration  of  vesicles  without  any  links  or  sandwiches 
attached . 

The  conservation  equations  can  be  derived  to  yield 
n  n-f 

V°  =  l  l  df  s(Kfi)f+s(Ksa)sV0  0  (17) 

f =0  s=0 


n 

Ab°  =  a  +  Ksal  +  X 

f=0 


n-f 

l  (s)df  s(Kfi )  f+s(Ksa)sV0  0 
s=0 

(18) 


n  n-f 

L°  =  i  +  KsaX  +  l  X  (f+s)df  s(Kf*)f+s(Ksa)sV0  0 

f=0  s=0 


(19) 


which  can  be  simplified  with  Equation  9  and  its  derivative  to 
the  formulas  below 


V°  =  (l+xy+y)nV0/0  (20) 

Ab°  =  a  +  KsaX  +  nxy(l+xy)n_1(l+y+xy)n_1 

( 1+y+xy-y/ ( 1+xy ) )V0f0 

(21) 

L°  =  A  +  KsaX  +  ny(l+xy)n_1(l+y+xy)n“l 

[ 1+x ( l+2y+xy-y/ ( 1+xy ) ]V0/0 


where  x  =  Ks3 

y  =  K  fA 


(22) 


A.  Sandwich  Model  Results 


The  sandwich  model  is  best  used  to  represent  possible 
applications  in  immunolysis  assays  [4].  A  vesicle  can  be 
made  with  antibody  receptors,  antigenic  links,  and 
antibody  ligands.  This  would  provide  more  flexibility  in 
the  types  of  systems  to  be  modelled.  In  Figures  9  and 
10  the  fraction  of  vesicles  with  at  least  one 
receptor/ link/ ligand  sandwich  is  shown  with  typical  values 
for  the  binding  constants  for  both  binding  steps,  Kf  and 
Ks.  The  value  of  8.3xl0“9  M  for  the  constants  L°  and  Ab° 
was  chosen  to  allow  for  comparison  with  previous  models. 

In  Waite  and  Chang 's  paper  on  the  monvalent  receptor  and 
monovalent  ligand  [14],  results  were  shown  for  varying 
values  of  n  (number  of  receptor  sites)  versus  fraction  of 
vesicles  with  one  or  more  antibodies  bound  to  the 
antigenic  receptors  (Figure  11) .  At  values  of  n  =  100  and 
Ab°  =  8.3xl0-9  M,  approximately  seven  percent  of  the 
vesicles  had  one  or  more  ligands  bound.  In  the  sandwich 
model,  with  the  link  concentration  held  at  the  value  of 
8 . 3xl0~9  M  it  takes  relatively  large  concentrations  of 
Ab°,  about  8.8xl0”6  M,  to  reach  the  same  seven  percent. 

The  sandwich  model  also  allows  the  antibody  concentration 
to  be  held  constant  and  the  link  concentration  to  be 
varied.  In  this  context  the  link  concentration  would  need 
to  be  about  7.9xl0-5  M,  with  the  [Ab°]  =  8.3xl0"9  M,  to 


Monovalent  receptorrMonovalent  ligand 
model.  Ab°=8 . 3xlC>  M.  Number  of  receptors 
(n)  vs.  fraction  with  antibody  ligands  attached 


Plot  from  Waite  and  Chang's  mono-mono  paper 
(used  by  permission) 


reach  seven  percent.  This  shows  that  the  sandwich  model 
is  less  sensitive  than  the  direct  binding  model.  This 
result  is  to  be  expected  since  this  model  requires  a  two 
step  process  to  be  completed  versus  a  one  step  process. 

The  plots  also  show  a  smooth  increase  in  the  number  of 
sites  "sandwiched"  up  to  a  maximum.  On  the  plot  which 
varies  Ab°,  the  maximum  is  one.  This  plot  will  not  drop 
back  off,  as  does  the  monovalent  receptor,  bivalent  ligand 
model,  shown  in  Figure  15  [15],  since  there  isn't  any 
competition  between  the  links  and  the  ligands  to  get  to 
the  receptor  cites.  There  is  no  competition  since  the 
links  are  specific  to  the  receptor  on  one  end  and  the 
ligand  on  the  other.  The  plot  of  varying  L° 
concentrations  reaches  a  maximum  at  approximately  eight 
percent.  This  is  logical  since  it  was  shown  that  the 
sandwiches  form  at  higher  concentrations  of  Ab°  on  the 
first  plot,  then  the  low  concentration  of  the  Ab°  in  this 
plot  lowers  the  possible  number  of  sandwiches.  The  number 
of  links  attached  does  not  affect  the  number  of  sandwiches 
after  a  point  if  there  are  not  enough  ligands  present  to 
form  the  second  step  of  the  process. 

The  values  of  105  for  the  binding  constants  are 
typical  [16]  but  can  vary  from  system  to  system.  This 
system  is  a  simple  two  step  process.  Changing  one  binding 
constant  by  some  orders  of  magnitude,  while  holding  the 
second  constant,  then  changing  the  second  constant 


while  holding  the  first  to  its  original  value  keeps  the 
results  constant,  as  is  shown  in  Figure  12.  It  is  not 
pictured  but  is  logical  that  an  increase  in  a  single 
binding  constant  or  increasing  both  would  increase  the 
relative  response. 

IV.  AGGLUTINATION  MODEL 


Another  model  which  can  be  ,sed  to  represent  the 
B-cell  response  is  the  monovalent  receptor,  multivalent 
antibody,  with  possibility  of  inter-linking  between 
vesicles,  called  agglutination  [4].  This  is  diagrammed  in 
Figure  13.  The  antibody  can  be  either  IgG  or  IgM.  The 
IgG  antibody  has  two  binding  sites,  while  the  IgM  has  ten 
binding  sites,  (1=2  for  IgG,  X  =  10  for  IgM,  where  X  is 
the  number  of  binding  sites  on  the  antibody)  but  for 
simplicity  in  this  model  it  will  be  assumed  that  only  two 
sites  will  be  available  for  binding  at  one  time.  Due  to 
steric  considerations  of  the  planar  shape  of  the  IgM,  this 
approximation  is  not  too  severe. 

From  Waite  and  Chang's  paper  on  multivalency  effects 
in  the  direct  binding  model  the  final  equilibrium  equation 
for  a  single  vesicle  with  monovalent  receptors  and 
multivalent  ligands  is 


Agglutination  Model 

Monovalent  receptor-bivalent  ligand 


A 


Figure  13 
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Vi,j  =  di/j(K1/a)j(XK1a)i+jv0/0  (23) 

where 

di,j  =  n!/[i! j 123 (n-i-2j) 1 ]  (24) 

and  V0  0  is  the  concentration  of  vesicles  without  any 
antibodies  bound,  either  as  danglers  (i  or  i')  or 
two-footers (j  or  j'),  Kj_  is  the  intrinsic  association 
constant,  K2  is  the  intra-cross-linking  constant,  and  a  is 
the  geometric  factor  which  takes  into  account  the  local 
concentration  enhancement  of  forming  a  two-footer  from  a 
dangler. 

Now  the  concentration  of  the  agglutinated  vesicles  can 
be  derived  from  this  information,  with  some  modifications. 
For  mathematical  simplicity,  the  number  of  vesicles  which 
can  form  inter-cellular  links  will  be  limited  to  two 
cells.  Therefore  the  degeneracy  of  the  agglutinated 
vesicles  can  be  shown  to  be  the  degeneracy  of  a  single 
vesicle  squared  with  the  number  of  bridges  (k)  taken  into 
account . 


dagglut 


n!n!/[i! j 1 2 3 (n-i-2j) !k!i' 


j  ' 123 ' (n-i'-2j  ') 1  ] 
(25) 


is  the  binding  constant  for  the  initial  bridge,  Kz 
is  the  binding  constant  for  the  bridges  following  the 
initial  bridge,  and  0  is  the  geometric  factor  which  takes 
the  steric  effects  and  local  concentration  effects 
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for  bridging  into  account.  Using  the  Law  of  Mass  Action 
as  has  been  shown  before,  the  equilibrium  equation  for  the 
concentration  of  two  agglutinated  vesicles  (C2)  can  be 
given  by 


c2  =  dagglut(K2/Q!)j(:vKla)i+jvolo(K2/a)j,(XKla)i,+j/vo,o 

(26) 


which  can  be  simplified  to 


c2  =  dagglutC^/^) j+j'(^Kla)i+i/+j+j/v2o,oKi(K2/^)k"1 

(27) 

Equations  23  and  27  can  be  combined  to  give  the  conservation 
equation,  which  has  been  abbreviated. 

(n)  n  (n-f) (n-f) 

V°  =  H  vi,  j  +2  l  H  H  C2  (28) 

i  j  k=l  i  j  i'j' 

Combining  similar  terms  and  substituting  allow  the 
simplified  equations  to  be  found. 

V°  =  G(n)V0f0  +  H(n)V20f0  (29) 


where 

n/2i 

G(n)  =  l  S(i)  (x1  +  xn-1)/(l+<Sn/2>i) 
i=0 


(30) 


The  n/2A  means  the  term  will  be  truncated,  and  the 
Kronecker  delta  function  <Sj^j  =  1  for  i  =  j  and  $i,j  -  0 
for  i»*j. 

i 

S ( i)  =  l  n!y!/[(i-j) ! j ! 23 (n-i-j ) !]  (31) 

j=0 

and 

n 

H(n)  =  X  (n!n!/[k! (n-k) ! (n-k) ! ] }xk2k-1Ki2G(n-k) 
k=l 

(32) 

x  =  aK^ 
y  =  K2/a 
z  =  Kz/ p 

The  antibody  conservation  equation  can  be  described  in 
similar  terms 


Ab°  =  a  +  J (n) VQ  0  +  W(n)V20/0  (33) 


J(n)  and  W(n)  are  defined  as 

n/2i  •  „  . 

J (n)  =  l  (ix1  +  (n-i)xn-1)S(i)/(l  +  «n/2,i) 
i=0 

(34) 


W(n)  =  l  (n!n!/k!  (n-k)  !  (n-k)  !  )xkzk-1Kj, 
k=l 

[kG2 (n-k)  +  2G(n-k) J(n-k) ] 

(35) 


and  x,  y,  and  z  are  defined  as  before. 


AGGLUTINATION  MODEL  RESULTS 


The  agglutination  model  is  useful  for  modeling  many 
systems.  It  can  be  used  to  model  the  immunolysis  assays 
with  monovalent  (Fab  fragment)  receptors  and  IgG  antibody 
or  bivalent  antigen  ligands  or  bridging  between  the  AIDS 
viruses  in  immunoadhesin  experiments  [8],  Figure  14  shows 
the  fraction  agglutinated  of  the  system  with  realistic 
binding  constants  for  the  formation  of  two-footers  (j)  and 
danglers  (i)  while  the  binding  constant  for  the  formation 
of  the  initial  bridge  is  varied.  This  is  done  since  no 
realistic  value  for  this  binding  constant  is  known.  The 
values  of  on  this  plot  vary  from  values  of  10”11  to 
103.  A  desired  value  would  be  to  find  the  antibody 
concentration  where  the  agglutination  would  be.  maximized. 
It  is  not  intuitive  that  the  maximum  be  consistent  for  all 
values  of  K^,  but  as  the  plot  shows,  no  matter  what  the 
probability  of  forming  the  first  bridge  is,  the  maximum 
occurs  at  the  same  concentration,  about  5.0xl0”6  M. 

From  Waite  and  Chang's  paper  on  multivalency  effects 
[15]  on  the  direct  binding  model  it  was  shown  that  the 
maximum  value  for  two-footers  is  about  5.0xl0-6  M  (shown 
in  Figure  15) .  This  maximum  is  the  same  for  fraction 
agglutinated,  which  implies  that  the  two  are  related.  In 
low  antibody  concentration  very  few  danglers  are  formed 
since  the  local  concentration  of  receptors  is  so  high  that 
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most  move  directly  to  two-footers.  A  bridge  cannot  be 
formed  if  both  sites  on  the  antibody  are  occupied, 
therefor  it  requires  danglers  to  be  present.  After  the 
maximum  concentration  of  two  footers  is  reached,  the 
relative  number  of  two-footers  per  vesicle  drop  due  to 
competition  with  danglers  for  receptor  sites.  These 
danglers  form  into  bridges.  The  fraction  agglutinated 
drops  off  when  the  antibody  concentration  gets  too  high. 
This  is  caused  by  all  of  the  vesicle's  receptors  being 
tied  up,  until  finally  there  are  no  open  receptors  on 
either  vesicle,  and  therefore  no  bridging. 

Figure  16  shows  the  effects  of  number  of  sites  per 
vesicle  (n)  upon  the  fraction  agglutinated.  This  curve  is 
relevant  since  the  computer  time  used  increases  rapidly 
with  increasing  n  and  realistic  numbers  of  receptor  sites 
would  take  too  long  to  calculate,  if  possible.  If  the 
effects  were  to  form  a  smooth  curve  the  effects  could  be 
extrapolated  to  a  realistic  number  of  receptors  per 
vesicle. 

This  plot  uses  an  antibody  concentration  of  10-12  M, 
to  set  a  low  value  to  try  and  approximate  the  threshold 
detection  value  for  different  methods.  The  light 
scattering  techniques  in  Cohen  and  Benedek's  papers 
[17,18]  rely  on  agglutination  for  detection.  Their  papers 
describe  laser  light  and  quasi-elastic  light  scattering 
from  the  early  stages  of  agglutination  for  low  level 
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Figure  16 


detection  of  ligands.  The  plot  could  also  be  related  to 
the  immunolysis  experiment  which  relies  on  the 
self -quenching  fluorescent  compound  carboxy fluorescein, 
and  the  lowest  level  of  light  detectable  by  the  light 
meter  to  show  a  response.  This  diagram  used  an  arbitrary 
fraction  of  agglutination  of  2xl0“9,  which  relates  to 
8.4xl0“22  M  of  the  vesicles  to  be  lysed,  to  estimate  the 
detection  threshold.  At  this  lower  limit  of  agglutination 
detection,  the  vesicles  would  need  at  least  675  receptor 
sites  to  detect  an  antibody  concentration  of  10”12  M. 


V.  NEW  DERIVATIONS 

Two  more  models  have  been  derived,  but  have  not  been 
programmed.  These  models  are  extensions  of  the  two  which 
have  been  presented  above. 

A.  Antibody  Multivalency  Effects  on  the  Sandwich  Model 

The  sandwich  model  can  be  made  more  realistic  by 
taking  into  account  antibody  multivalency.  As  before,  the 
IgG  will  be  represented  ideally  and  the  IgM  will  be 
limited  to  an  effective  binding  number  of  two  to  limit  the 
mathematical  complexity. 


This  model  again  uses  the  vesicle  with  n  receptor 
sites  which  can  be  involved  first  with  the  linking  ligand, 
then  the  antibody.  The  degeneracy  of  this  system  is 

df,s,t  =  n!/[f !s!t!2t(n-f-s-2t) ! ]  (36) 

Figure  17  shows  this  model  where  f  receptor  sites  have 
a  link  attached,  s  with  an  antibody  bound  to  a  link 
singularly,  in  a  similar  manner  as  a  dangler,  and  t  with 
links  bound  by  the  same  antibody. 

The  Law  of  Mass  Action  equations  necessary  to  derive 
the  equilibrium  equations  are 

$ 

Kt  t  times 

v/f,s,t  *  v'f,s+l,t-l  •••  vf,s+t,o  (3?) 


Ks  s+t 

v#f,s+t,o  “  v'f+l,s-l,o  +  a'  vf+s+t,o,o  +  a 

times 

(38) 


Kf 


f+s+t 


f+s-ft,o,o 


v#f+s+t-l,o,o  +  * ' 


+  J i 


times 


(39) 
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These  equations  taken  to  the  limits  shown  above  yield 


the  equilibrium  equation 

vf,s,t  “  <3f,s,t(Kt/»)t(Ksal)s+t(Kfl)f+s+tv0_0<0 

(40) 

1  is  the  number  of  receptor  sites  on  each  antibody  and  a 
is  the  geometric  factor  which  takes  local  concentration 
effects  into  account. 

Defining  x  =  Kfi ,  y  =  Ksa  ,  and  z  =  Kt/a  the 
following  conservation  equations  can  be  derived. 


n-f-s/2 
V°  —  l 
t=0 


nif  Z  df ,  s ,  txf +s+tys+t z ^o  ,0,0 

S=0  f =0 


(41) 


n-f-s/2  n-f 
Ab °  =  I  l 

t=0  S=0 


n 

l 

f=0 


z  df,s,t(s+t)xf+s+tys+tztvo,o,o 


(42) 


n-f-s/2  n-f 

L°  =  l  l 

t=0  S=0 


Z  df , s , t ( f +S+2 t ) xf +S+tyS+tztVo 


Using  Equation  9  and  its  derivative  the  conservation 


equations  can  be  simplified  into 


V°  =  ( 1+x) n( 1+xy) n-fH(t) 


(44) 


Ab°  =  a  +  2Ksai  +  K2sa2i  + 

(n-f ) xy ( 1+x) n ( 1+xy) ( t) v0f  0  f Q 

(l+x)n(l+xy)n-ftH(t)V0>0f0 

(45) 

L°  =  i  +  2Ksai  +  K2sa2X  +  nx(l+x)n“1(l+x)n~fH(t)V000 

+  ( 1+x) n (n-f) xy ( 1+xy) n“f_1H (t) VQ  f0/0 
+  (l+x)n(i+xy)n-f2tH(t)V0/0/0 

(46) 

where 

n-f-s/2 

H (t)  =  l  (n-f-s) !/[t!2t(n-f-s-2t) !]  (xyz)t 
t=0 

(47) 

The  results  from  this  model  would  be  expected  to  be 
related  to  the  monovalent  ligand  sandwich  model  much  in 
the  same  manner  as  Waite  and  Chang's  papers  on  monovalent 
direct  binding  and  multivalency  effects  are  related 
[14,15].  It  would  be  expected  that  these  effects  include 
geometric  factors  for  local  concentration  of  receptors 
after  one  end  of  the  antibody  is  attached.  It  is  also 
logical  that  the  number  of  sandwiches  with  both  sites  on 
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the  ligand  being  attached  having  a  peak  when  compared  to 
initial  antibody  concentration.  This  would  be  caused  by 
the  competition  of  the  "two-footed  sandwiches”  with  the 
danglers  for  receptor  sites  with  links  attached. 

B.  Bivalent  Receptor/Bivalent  Ligand  Model 

The  bivalent  receptor,  bivalent  ligand  model  is  the 
most  complex  model  presented.  The  vesicle  surface  is 
covered  with  n  surface  immunoglobulins  (slg)  which  each 
have  two  receptor  cites.  This  system  is  shown  in  Figure 
9. 

Theoretically  a  B-cell  has  antibody  chains  from 
surface  Ig  to  slg  which  are  up  to  twenty  receptors  long. 
These  chains  can  form  loops  upon  themselves,  long 
continuous  chains,  or  large  clusters.  The  number  of 
possible  arrangements  is  mathematically  intractable.  To 
make  this  model's  derivation  somewhat  straightforward  this 
model  limits  the  chaining  effects  to  two  slg.  The 
possible  arrangements  of  two  receptor  cites  is  also  shown 
on  Figure  9. 

The  degeneracy  of  this  system  need  to  take  this  2n 
factor  into  account,  but  this  falls  out  since  all  possible 
combinations  of  antibiotic  attachment  and  empties  also 
have  this  2X  factor. 
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da_g  =  n!/[a!b!c!d!e! f !g! (n-a-b-c-2d-2e-2f-2g) ! ] 

(48) 


The  Law  of  Mass  action  equations  are  shown 
diagrammatical ly  on  Figures  19a  and  19b. 

The  breakdown  of  each  structure  is  not  always  a 
singular  path,  but  this  system  falls  under  the  same  rules 
as  a  state  function,  thus  any  one  pathway  may  be  selected 
for  each  combination.  The  one  pathway  chosen  must  be 
consistent  for  all  the  Mass  Action  equations.  The 
equilibrium  equation  for  any  specific  combination  of 
antibodies  is  given  by 


va-g  =  da_g(Kaab)a+b+c+d+2e+2f+29(Kbab)b(Kc)c(Kd)d 

(Ke)e+f+g(Kf)f(Kgab)9 


(49) 


This  equation  uses  ab  to  represent  free  antibody 
concentration  to  eliminate  confusion  between  it  and 
receptor-ligand  combination  a.  The  binding  constants  Kb 
through  Kg  all  have  geometric  local  concentration  effects 
tied  into  their  values. 
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The  conservation  equation  for  the  vesicles  is 


£  <p  e  u  f  n-a  n 

v°  =1111111  ^a-g  (Kaab) a  (KaKbab2 ) b 
g=0  f=0  e=0  d=0  c=0  b=0  a=0 

(KaKcab)  c(KaKdab)  d  (K2  aKeab2  )  e  (K2  aKeKf ab2  )  f  (K2  aKeKgab3 )  9 

f  =  n-a-b 
u  =  n-a-b-c/2 
e  =  n-a-b-c-2d/2 
<t>  =  n-a-b-c-2d-2e/2 
£  =  n-a-b-c-2d-2e-2f/2 

(50) 


Models  have  been  derived  with  similar  arrangements 
(21]  of  bivalent  ligands  and  bivalent  receptors,  but 
again,  this  model  is  receptor-carrier  based  [10].  This 
will  allow  questions  about  linked  receptors  per  vesicle  to 


be  answered. 


VI.  CONCLUSION 


The  purpose  of  this  paper  was  to  describe  in 
detail  the  derivation  of  a  theoretical  model  for 
antigen-antibody  interaction  and  how  that 
interaction  relates  to  the  immune  response, 
immunoassay  methods  of  detection,  and  the  use  of 
immunoadhesins  for  AIDS  therapy.  These  models  were 
receptor  carrier  based  which  allows  the  questions 
of  number  of  interactions  per  vesicle  or  cell  to  be 
addressed.  The  results  may  be  used  to  discover 
parameters  for  these  systems.  With  help  from 
experimental  immunologists  to  relate  real  known 
binding  constants  and  concentrations,  boundary 
limits  and  interactions  on  the  cell/vesicle  surface 
may  be  described.  These  models  can  be  used  as 
building  blocks  for  describing  more  complex 
systems . 
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APPENDIX  1 


Computer  programming  of  these  equations  does 
pose  some  interesting  problems.  The  precision 
for  the  VAX  is  normally  e~88  and  double  precision 
is  e-300.  One  of  the  problems  occurs  when  adding 
a  series  of  numbers  of  greatly  different  orders 
of  magnitude.  Even  when  working  in  natural 
logarithims,  the  numbers  must  be  changed  back  to 
their  non  logarithmic  form  to  add  them.  This 
problem  can  be  avoided  by  using  a  series  of  steps 
to  make  the  numbers  more  manageble. 

Desired  quantity  In  (x3  +  x2  +  x3) 


(a) 

Find:  In  x1#  In  x2,  In  x3 

(b) 

Find  largest  quantity  from  (a) 

(c) 

Subtract  (b)  from  all  values  in 

(a) 

(d) 

Take  the  exponential  of  values 

in  (c) 

(e) 

Sum  numbers  in  d  (magnitudes  close) 

(f) 

Take  In  (e) 

(g) 

Add  In  (f)  +  In  (b) 

(h) 

Equals  In  (X!  +  x2  +  x3) 

i .  e. 

Desired  104  +  10-1  (Log10  will 

be  used) 

(a) 

log  104  =  4,  log  10-1  =  -1 

(b) 

largest  =  4 

(c) 

0,-5 

(d) 

10°  =  1,  10“5  =  .00001 
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(e)  1.00001 

(f)  log  1.00001  -  4.3X10"6 

(g)  4  +  4.3X10-6  =  4.0000043 

(h)  lo(9)  =  10000.1  as  expected 

The  other  complication  arises  when  the 
factorial  for  large  numbers  is  needed.  This 
arises  in  all  the  degeneracy  calculations  since 
n!  is  a  very  large  number.  This  can  be 
approximated  using  Stirling's  approximation  which 
says: 

ln(n! )  =  (n+1-1/2) ln(n+l)  -  (n+1)  + 
l/21n(2?t)  + 

1/12 (n+1)  -  1/360 (n+1) 3  +  1/1260 (n+1) 5  ... 

This  approximation  gives  very  exact  answers 
without  overloads  '•'he  computer.  Most  numbers 
would  overflow  normal  precision  since  100!  = 
9.3xl0157,  and  many  would  averflow  double 
precision,  therefore  Stirling's  approximation  is 
used  whenever  n  is  greater  than  seven. 


Appendix  2 

This  computer  program  is  used  to  develop  the  plots  and 
results  for  the  sandwich  model.  The  loops  fill  a  two- 
dimensional  array  for  all  VfiS  combinations.  This  array  is 
forced  into  a  one-dimensional  array  where  the  summation  of 
the  fraction  can  be  found.  The  data  is  then  passed  to  data 
files  to  be  interpreted. 
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ccc  PROGRAM:  SAND. FOR 

c  To  run  type: 

c  fortran/g_f loating  sand. for 

c  link  sand 

c  run  sand 

c 

implicit  real*8  (a-h,o-z) 

dimension  Vfs(200,200) 

dimension  temp (40000) 

dimension  temp2 (40000) 

akf=l.d5 

aks=l.d5 

ictr=l 

nsand=l 

sum=0 

Vo=dlog ( 4 . 2d-13 ) 
do  5  n=100, 101, 10 
do  10  ih=900 , 901 , 2 
a=8 . 3*10. ** (-ih/ 100 . ) 
do  15  im=900, 300,-50 
al=10. ** (-im/100 . ) 
x  =  aks*a 
y  =  akf*al 
w  =  aks*a*al 
Voo  =  Vo-n*dlog(l+x*y+y) 
do  20  if  =  0,n 
do  25  is  =  o,n-if 

Dfs=fact (n) -fact (if ) -fact (is) -fact (n-is-if ) 
Vfs(if+1, is+1) =Dfs+ (if+is) *dlog(y) +is*dlog (x) +Voo 
25  continue 
20  continue 
nn=n-l 
yy=l+x*y 

aLo=Voo+dlog (n*y ) +nn*dlog (yy ) 


alLo=Voo+dlog (n*x*y) +nn*dlog (y+yy) 
a2Lo=dexp(aLo) +al+w+dexp (alLo) 
goto  ll 

Abo=dlog (n*x*y ) +nn*dlog (y+yy) +dlog ( 1-y /yy ) +V00 
Abo2=a+w+dexp (Abo) 

11  continue 
ictr=l 

do  100  ia=0,n 
do  95  ib=0,n-ia 
if  (ib. lt.nsand)  goto  95 
Temp ( ictr ) =Vf s ( ia+1 , ib+1 ) 
ictr=ictr+l 
95  continue 
100  continue 

do  110  ic=l, (ictr-1) 
big=temp(l) 

if  (temp ( ic) .gt. big)  big=temp(ic) 

110  continue 
sum=0 

do  120  id=l, (ictr-1) 
temp2(id)  =  temp (id) -big 
temp2 (id)=dexp(temp2 (id) ) 
sum=temp2 ( id) +sum 
120  continue 

sum=sum*dexp (big) 
write(25, 999)  n 

999  format ('  current  n  value  is:  ',i5) 

write(25, 1000)  Abo2 

1000  format ('  computed  value  of  Abo  is:  ',dl6.8) 

write (25, 1001)  a2Lo 

1001  format {•  • computed  value  of  Lo  is:  ',dl6.9) 
f raction=sum/dexp (vo) 

write(25, 1002)  nsand, fraction 

1002  format ('  ' fraction  with 
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l',i5,'or  more  sandwiches  is:  ',d!6.8) 
write(20, 1003)  dloglO(a2Lo) , fraction, dloglO (fraction) 
write (21, 1004)  n, fraction, dloglO ( fraction) 

1003  format (3el6. 8) 

1004  format(i5,2el6.8) 

15  continue 

10  continue 
5  continue 
stop 
end 

real*8  function  fact(n) 
implicit  real*8  (a-h,o-z) 
tl=. 5*dlog (2 .d0*3 . 14 1592 653 58 97 9d0) 
if  (n.gt.6)  goto  30 
if  (n.eq.0)  fact=0.0 
if  (n.eq.l)  fact=0.0 
if  (n.eq.2)  fact*dlog(2.d0) 
if  (n.eq.3)  fact=dlog(6.d0) 
if  (n.eq.4)  fact=dlog(24 .d0) 
if  (n.eq.5)  fact=dlog(120.d0) 
if  (n.eq.6)  fact=dlog(720.d0) 
return 
30  z=n+l 

fact=(z-. 5) *dlog (z)-z+tl+l./(12.*z)-l./(360. 
1*Z**3)+1./ (1260.* 

1Z**5)-1./(1680.*Z**7) 

return 

end 
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Appendix  3 

This  program  is  more  complex  and  layered  to  solve  the 
agglutination  model.  The  functions  S(i),  G(n)  and  others 
which  are  used  more  than  once  have  their  own  arrays  to  make 
the  calculations  easier. 
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ccc 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


PROGRAM:  BEST. FOR 

to  run,  type 

fortran/g_f loating  [waite. inter ] best. for 

link  best 

run  best 

type  forOOO.dat 


Philosophy  of  program: 

sums  of  large  numbers  are  handled  with  the 
logarithmic,  reduce  by  biggest,  exp,  sum, 
log,  increase  by  biggest  method, 
the  quadratic  formula  is  implemented  with  similar 
consideration  given  to  large/ small  numbers  in 
function  f 

Sterling's  formula  is  implemented  in  function  fact 


The  summary  of  variables  possessing  the  various  functions  is  as 
follows: 


(nsize) 

s (nsize+l, k+1)  is  the  s(i) 
g(nsize+l)  is  the  g(nsize) 
j (nsize+l)  is  the  j (nsize) 
hn  is  the  h(n) 
wn  is  the  w(n) 

pn  is  involved  with  the  m  or  more  bridges  calculation 


implicit  real*8  (a-h,o-z) 

dimension  tl(300) ,bigest(300) ,suml(300) , s (300, 300) 
dimension  zz(300,300) , zzn ( 300 , 300) ,sum(300) ,tut(300) 
dimension  uu (300, 300) , uun(300, 300) , sumr  100) ,sl(300) 
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dimension  s2  (300) , s3  (300) 

zero=0 . 0 

mbrdg=3 

do  987  qw— 5,7, 1 
akl=10. ** (qw) 
ak2=l.d5 
aki=l.d3 
akz=l . d5 
val=2 . 0 
alpha=7 .d5 
beta=7 .d5 
y=ak2 /alpha 
z=akz/beta 
vo=4.2d-13 
write (0,2323) 

2323  format ( '  ' , '  kl  k2  ki  kz 

write (0,1212)  akl , ak2 , aki , akz , val 

1212  format ('  ',5dl6.8) 
write (0, 3434) 

3434  format ('  alpha  beta  vo') 

write (0,1212)  alpha,  beta,vo 
do  7860  n=100, 101, 10 
write(0,8642) 

8642  format ('  ' ) 

write(0,4545)  n 

4545  format ('  'current  value  for  n  is  :  ',i6) 
do  5  nsize=0,n 
do  15  k=0,nsize/2 
do  25  j=0,k 

tl ( j+1) =f act (nsize) -fact (k-j ) 

1-f act ( j ) -f act (nsize-k- j ) +j*dlog (y/2 . ) 

25  continue 

bigest (k+1) =tl (1) 
do  35  j=0,k 


val' ) 


if (tl ( j+1) .gt .bigest (k+1) )  bigest (k+1) =tl ( j+1) 

35  continue 

suml (k+1) =0 . 0 
do  45  j=0,k 

tl ( j+1) =tl ( j+1) -bigest (k+1) 
if (tl ( j+1) . It . -300 . )  goto  102 
tl( j+l)=dexp(tl( j+1) ) 
goto  101 
102  tl(j+l)=0.0 

101  suml  (k+1) =suml (k+1) +tl (j+1) 

45  continue 

s (nsize+1 , k+1) =dlog (suml (k+1) ) +bigest (k+1) 

15  continue 
5  continue 

do  666  i=800, 300,-5 
a=10 . ** (-i/100 . ) 
write(0, 8642) 
write(0, 5676)  a 

5676  format ('  ' , '  value  for  free  antibody  a  is  :',dl6. 
x=a*val*akl 
do  135  nsize=0,n 
do  145  k=0,nsize/2 
tll=k*dlog(x) 
t2=(nsize-k) *dlog(x) 
if(k.eq.O)  t3=-10000. 
if(k.eq.O)  goto  17 
ak=k 

t3=tll+dlog (ak) 

17  if (nsize . eq.k)  t4=-10000. 
if (nsize. eq.k)  goto  18 
anmk=nsize-k 
t4=t2+dlog (anmk) 

18  if ( (nsize/2) *2 . eq. nsize)  goto  3 
goto  333 


3  if (k.eq. (nsize/2) )  tll=tll-dlog (2 .do) 
if (k. eq. (nsize/2) )  t2=t2-dlog (2 .do) 
if (k.eq. (nsize/2) )  t3=t3-dlog(2 .do) 
if (k.eq. (nsize/2) )  t4=t4-dlog(2 .do) 

333  zz (nsize+1, k+1) =s (nsize+1, k+1) +tll 
zzn(nsize+l,k+l) =s( nsize+1, k+1) +t2 
uu(nsize+l,k+l)=s ( nsize+1, k+1) +t 3 
uun(nsize+l,k+l)=s(nsize+l,k+l)  +  t4 
145  continue 

big=zz (nsize+1, l) 
bigu=uu (nsize+1 , 1) 
do  56  k=0, nsize/2 

if (zz (nsize+1, k+1) .gt.big)  big=zz (nsize+1, k+1) 
if (zzn(nsize+l,k+l) .gt.big)  big=zzn (nsize+1 , k+1) 
if (uu(nsize+l,k+l) .gt.bigu)  bigu=uu (nsize+1 , k+1) 
if (uun(nsize+l,k+l) .gt.bigu)  bigu=uun (nsize+1 , k+1) 
56  continue 

do  67  k=0, nsize/2 

zz (nsize+1, k+1) =zz (nsize+1, k+1) -big 
zzn( nsize+1, k+1) =zzn( nsize+1, k+1) -big 
uu( nsize+1, k+1) =uu ( nsize+1, k+1) -bigu 
uun (nsize+1, k+1) =uun (nsize+l , k+1) -bigu 
67  continue 

sum (nsize+1) =0 . 0 
sumu (nsize+1) =0 . 0 
do  78  k=0, nsize/2 

if (zz (nsize+1, k+1) . It. -300. )  goto  578 
zz (nsize+1 , k+1) =dexp ( zz (nsize+1 , k+1) ) 
goto  478 

578  zz (nsize+1 , k+1) =0 . 0 
478  continue 

if (zzn(nsize+l,k+l) . It. -300. )  goto  579 
zzn( nsize+1, k+1) =dexp ( zzn (nsize+1 , k+1) ) 
goto  479 
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479  continue 

if (uu(nsize+l,k+l) .It. -300. )  goto  679 
uu (nsize+1, k+1) =dexp(uu (nsize+1, k+1) ) 
goto  680 

679  uu(nsize+l,k+l)=0.0 

680  continue 

if (uun(nsize+l,k+l) .It. -300. )  goto  779 
uun (nsize+1, k+1) =dexp (uun(nsize+l,k+l) ) 
goto  780 

779  uun (nsize+1, k+1) =0.0 

780  continue 

sum(nsize+i) =sum (nsize+1) +zz (nsize+1, k+1) 
l+zzn( nsize+1, k+1) 

sumu (nsize+1) =sumu(nsize+l) +uu (nsize+1, k+1) 

1+uun (nsize+1 , k+1) 

78  continue 

sum (nsize+1 )=dlog (sum (nsize+1) )+big 
sumu(nsize+l)=dlog(sumu(nsize+l) ) +bigu 
135  continue 

do  601  k=l,n 

si (k) =2 . *fact (n) +k*diog(x) + (k-1) *dlog (z) +dlog (aki) -fact (k) - 
1  2.*fact(n-k) 
ak=k 

s2 (k)=sl(k)+dlog(ak)+2 . *sum(n-k+l) 
s3 (k)=sl(k)+dlog(2 .d0)+sum(n-k+l)+sumu(n-k+l) 
si (k) =sl (k) +dlog(2 . dO) +2 . *sum(n-k+l) 

601  continue 
bl=sl (1) 
b2=s2(l) 

do  602  k=l,n 
if (si (k) .gt.bl)  bl=sl (k) 
if (s2 (k) .gt .b2)  b2=s2 (k) 
if (s3 (k) .gt .b2)  b2=s3(k) 

602  continue 
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si (k) =sl (k) -bl 
s2  (k)  =s2  (k)  -b2 
s3(k)=s3(k)-b2 

603  continue 
hn=0 . OdO 
wn=0. OdO 
pn=0. OdO 

do  604  k=l,n 

if (sl(k) .It. -300. )goto  720 
sl(k)=dexp(sl(k) ) 
goto  721 

720  si (k) =0 . 0 

721  continue 

if (s2 (k) . lt.-300. ) goto  722 
s2 (k)=dexp(s2 (k) ) 
goto  723 

722  s2 (k) =0 . 0 

723  continue 

if (s3(k) .It. -300.)  goto  724 
s3 (k)=dexp(s3 (k) ) 
goto  725 

724  S3 (k) =0 . 0 

725  continue 
hn=hn+sl (k) 
wn=wn+s2 (k) +s3 (k) 

if (k. lt.mbrdg)  pn=pn+sl(k) 

604  continue 
pn=hn-pn 
hn=dlog(hn) +bl 
wn=dlog(wn)  +  b2 
pn=dlog(pn) +bl 
bbb=sum(n+l) 
ccc=dlog (vo) 
aaa=hn 


bog=aaa 

i f ( bbb . gt . bog ) bog=bbb 

if ( ccc . gt . bog) bog=ccc 

aaaa=aaa-bog 

bbbb=bbb-bog 

cccc=ccc-bog 

vOO=f (aaaa, bbbb, cccc) 

8889  if (vOO .gt. -300 .) write (0, llll)  dexp(vOO) 
if (v00.le.-300. ) write (0, llll)  zero 
terml=dlog(a) 
term2=sumu (n+1) +v00 
term3=wn+2 . *v00 
big=terml 

if (term2 .gt.big)  big=term2 
if (term3 .gt.big)  big=term3 
terml=terml“big 
term2=term2-big 
term3=term3  -big 
stum=0 . 0 

if (terml. It . -300 . 0) goto  851 
termi=dexp ( terml ) 
goto  852 

851  terml=0.0 

852  continue 

if (term2 . It. -300. )  goto  853 
term2=dexp ( term2 ) 
goto  854 

853  term2=0.0 

854  continue 

if (term3 . It. -300. )  goto  855 
term3  =dexp ( term3 ) 
goto  856 

855  term3=0.0 

856  continue 


856  continue 

abo=terml+term2+term3 
abo=dlog ( abo) +big 

if (abo.gt. -300. )  write(0, 1113)  dexp(abo) 
if (abo. le. -300. )  write (0, 1113)  2ero 
write(0, 1114)  dloglO (dexp (abo) ) 
guppy=hn+2 . *v00-ccc 

if (guppy .gt. -300. )  write(0,222)  dexp(guppy) 
if (guppy . le. -300. )  write(0,222)  zero 
patch=pn+2 . *v00-ccc 

if (patch. gt. -300. )  write(0,242)  mbrdg, dexp (patch) 
if (patch. le. -300. )  write(0,242)  mbrdg, zero 
write ( 1 , 111)  dloglO (dexp (abo) ) , dexp (guppy) 
write (2, 111)  dexp (abo) , dexp (guppy) 
write(3, 111)  dloglO (dexp (abo) ) ,dexp(patch) 
write (4 , 1121)  n, dexp (guppy) 

242  format  ( '  ' fract.  with  ',i5,'  or  more 
lbridges  :  ',dl6.8) 

111  formats  ',3el6.8) 

1121  format ( '  ',i5,3el6.8) 

1111  format ('  ' computed  value  of  vOO  is:  ',dl6.8) 

1113  format ('  computed  value  of  abo  is:  ',dl6.8) 

1114  format ( '  ' computed  value  of  log (abo)  is:  ',dl6. 

222  format^  ' f(agglut)=  ',dl6.8) 

666  continue 
7860  continue 
987  continue 
stop 
end 

real*8  function  fact(n) 
implicit  real*8  (a-h,o-z) 
tl=.5*dlog(2.d0*3. 14159265358979d0) 
if(n.gt.6)  goto  10 
if(n.eq.O)  fact=0.0 


if(n.eq.l)  fact=0.0 
if(n.eq.2)  fact=dlog(2 .dO) 
if(n.eq.3)  fact=dlog(6.d0) 
if(n.eq.4)  fact=dlog(24 .dO) 
if(n.eq.5)  fact=dlog(120.d0^ 
if(n.eq.6)  fact=dlog(720.d0) 
return 
10  z=n+l 

fact= ( z- . 5) *dlog(z) -z+tl+1./ (12 . *z) -1./ 
l(360.*z*z*z)+l./(1260.* 
lz*z*z*z*z)  -l./(i680.*z*z*z*z*z*z*z) 
return 
end 

real*8  function  f(a,b,c) 
implicit  real*8  (a-h,o-z) 
x=dlog ( 4 . do ) +a+c-2 . *b 
111  format ('  ',3dl6.8) 
if(x.gt.O.O)  goto  15 
if (x.lt.-5.0)goto  5 
x=dexp(x)+l. 

x=dsqrt ( x ) *dexp ( b ) -dexp ( b ) 
f=dlog(x) -dlog(2 .dO) -a 
return 

5  tl=dexp(a+c-b-b) 

f=C-b-tl+2 . *tl*tl-5 . *tl*tl*tl+14 . *tl**4-42 . *tl* 
return 
15  y=2*b-a-c 

if (y. It. -300. )goto  35 

y=dexp(y) 

goto  36 

35  y=0 . 0 

36  continue 

z=(-dsqrt (y)+dsqrt (y+4 .dO) )/2. 
f=dlog ( z ) + . 5*c- . 5*a 


